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1. Introduction 



The study of heavy-ion collisions at intermediate energies (50 < E < 
1000 MeV/nucleon) provides a rich source of information for many rare phenomena 
QQ ■ such as multifragmentation, collective flow as well as particle production [1,2]. One can 

also shed light on the mechanism behind the fragmentation in highly excited nuclear sys- 
CO ■ tems. In this energy region, multifragmentation appears to be a dominant de-excitation 

lO I channel apart from the other less populated channels of manifestation of liquid gas phase 

ly^ . transition in finite nuclear systems [1,3,4]. In the literature, multifragmentation has also 

C^ ' been considered as a gateway to nuclear equation of state [5,6]. Numerous investigations 

^D , are cited in the literature which handle the de-excitation of nuclear system in multifrag- 

mentation [7-13]. 

The experimental analysis of the emission of intermediate mass fragments (IMF's), (5 < 
, , A < Atot/^), has yielded several interesting observations: De Souza et al. [10] observed 

r> \ a linear increase in the multiplicity of IMF's with incident energies for central collisions. 

t^ . In this study, incident energy was varied between 35 and 110 MeV/nucleon. On the other 

hand, Tsang et al. [II] reported a rise and fall in the production of IMF's. The maximal 
value of the IMF's shifts from nearly central to peripheral collisions with the increase in 
the incident energy. More refined results were reported by Peaslee et al. [12] for the reac- 
tion of ae-RTr*'* + jgAw^^^ for incident energy between 35 and 400 MeV/nucleon. Their 
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analysis revealed that IMF's multiplicity first increases till 100 MeV/nucleon and then de- 
creases slowly. These findings pose a stringent test for any theoretical model designed for 
the study of multifragmentation. 

Theoretically, multifragmentation can be studied by statistical [2] as well as semi- 
classical [5] models, respectively. The relation between multifragmentation process and 
nuclear equation of state was extensively studied by several authors with in the statistical 
approach for intermediate energy heavy-ion collisions [2,4,6]. On the other hand, semi- 
classical dynamical models [5] are very useful for studying the reaction from the start to 
final state where matter is fragmented and cold. In addition, these models also give pos- 
sibility to extract the information about the nuclear equation of state [13] and NN cross- 
section [1,5,14]. The interaction among nucleons (in a heavy-ion reaction) can be studied 
within the G-matrix, with its real part representing the mean field and complex part denotes 
the NN cross-section [1,5,14]. Note that the contribution of imaginary part of the interac- 
tion is nearly absent in the low energy process such as fusion, fission and radioactivity 
[15]. One often uses a parametrized form for the real and imaginary parts of the G-matrix. 
It is well accepted to use a density dependent Skyrme type interactions for the real part 
of the G-matrix. However, heavy-ion dynamics depends not only on the density but also 
on the entire momentum plane [7]. Therefore, it is advisable to use momentum dependent 
interactions additionally. 

The exact nature of NN cross-section, on the other hand, is still an open question [7,14]. A 
large number of calculations exist in the literature suggesting different strength and forms 
of the NN cross-sections [14,16]. In a simple assumption of hard core radius of NN poten- 
tial, one has often used a constant and isotropic cross-section of 40 mb. In other calcula- 
tions, a constant and isotropic cross-section with magnitude between 20 and 55 mb is also 
used [14]. The most sensitive observable to pin down the NN cross-section is collective 
flow. Recent calculations advocated its strength between 35-40 mb [14]. We shall con- 
centrate here on multifragmentation. Our present study will be based on the semi-classical 
model, namely, quantum molecular dynamics (QMD). In semi-classical model, one has to 
face the problem of the stability of nuclei. Several model ingredients such as width of the 
Gaussian (L), has been used as free parameter. It varies between 4.33 fm? to 8.66 frn? 
[5]. On the same time, one needs to identify the clusters with the help of clusterization 
algorithm. Different cluster recognisation algorithms can also influence the fragmentation. 
Our present interest is to perform a comparative study of different model ingredients alogn- 
with different NN cross-sections and to see whether it is possible to pin down the strength 
of NN cross-sections by QMD model or not. 

Our article is organized as follow: we discuss the model briefly in section-II. Our results 
and discussions are given in section-Ill and we summarize the results in section-IV. 



2. Quantum Molecular Dynamics (QMD) Model 

The QMD model [1,5,7,14,17] is a time dependent N-body theory which simulates the time 
evolution of heavy-ion reactions on an event-by-event basis. It is based on the generalized 
variational principle. As with every variational approach, it requires the choice of a test 
wave function $. In the QMD approach, this is an N-body wave function with 6N time- 
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dependent parameters if nuclear system contains N nucleons. 

The basic assumption of QMD model is that each nucleon is represented by coherent states 
of the form (we set 7i, c = 1) which are characterized by 6 time dependent parameters, t\ 
and Pi, respectively. 

/ 2 \3/4 

Mr,t)=(—j cxp-(^~-'"')'/^cxp'(^"P"''W), (1) 

Here L is the width of the Gaussian distribution. This width varies between 4.33 and 8.66 
fm^ in the literature [5]. The total N-body wave function is assumed to be a direct product 
of the coherent states 

$ = nfj+-^-0,(f,rl,pl,t) • (2) 

To calculate the time evolution of a system, we start out from the action 

S^ f £[<P,$*]dt (3) 



with the Lagrange functional 

£ = {<P\ifi^^-H\<P) ■ (4) 

The total time derivative includes the derivation with respect to parameters. The time 
evolution of parameters is obtained by the requirement that the action is stationary under 
the allowed variation of wave function. This leads to a Euler-Lagrange equation for each 
time-dependent parameter. 

Thus, the variational principle reduces the time evolution of N-body Schrodinger equation 
to the time evolution equations of Q{At + Ap) parameters to which a physical meaning 
can be attributed. The equations of motion for the parameters pi and fl reads as: 

dfl _ d < H > dpi _ d < H > 

dt dpi ' dt dfl 

If < H > has no explicit time dependence, QMD conserves the energy and momentum by 

construction. 

The nuclear dynamics of the QMD model can also be translated into a semi classical 

scheme. The Wigner distribution function fi of a nucleon i*^ can be easily derived from 

the test wave functions. 

af,p,t) = _i_cxp-(^^-'"'(*))'-2/^.cxp-(p-P'(*))'-^/2''' (6) 

and the total Wigner density is the sum of those of all nucleons. Hence the expectation 
value of total Hamiltonian reads 

{H)={T)+{V) 

i i j>i 

xfj{f,p, t)dfdr' dpdp' (7) 
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Where V^^ is given as: 

^ V t' ) — ^Skyrme ^ '^Yukawa ^ ^Coul ^ ^ mdi 

h5{^-r)+h6{7-r)p^-^ U^ 

exp(-\r-r\/fi) ZjZjG^ 
(|r' — r|//i) |r — r| 

-^U\ll\tS^-vf + \W-T^ (8) 

which consists of Skyrme, Yukawa, Coulomb and momentum dependent parts of the inter- 
action and Zi, Zj are the charges of i*'' and j*'' baryons. 

Further, potential part resulting from the convolution of the distribution function /^ and ^ 
with the Skyrme interactions Vskyrme reads as : 

Vskyrn^e - "^ (^) + ^ (^Y " ^9) 

The two of the three parameters of equation of state are determined by demanding that 
at normal nuclear matter density, the binding energy should be equal to 16 MeV. The 
third parameter 7 is usually treated as a free parameter. Its value is given in term of the 
compressibility: 

2 92 fE\ 
- = 9pV [a) ■ ^''^ 

The different values of compressibility give rise to Soft and Hard equations of state. 
The momentum dependence (Vmdi) of NN interaction, which may optionally be used in 
the model, is fitted to experimental data [5] on the real part of the nucleon optical potential. 
This yields: 



Vm 



6.ln^e-{Ap)' + l)-(^ • (11) 



The potential part of the equation of state resulting from the convolution of distribution 
functions fi and fj with Skyrme and momentum dependent interactions, reads as: 

\ Po J \ Po 

+6-ln^e-{Ap)^ + l)-f^ ■ (12) 

For the value of constants and details, reader is referred to ref. [5]. The inclusion of MDI 
is labeled as Soft momentum dependent (SMD) and Hard momentum dependent (HMD) 
equations of state. In recent studies, momentum dependent forces are found to be essential 
part of dynamics. 
During the propagation, two nucleon can collide if they satisfy the condition \fi — fjl < 
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\J a JTx. In this procedure, a binary collision is blocked if the final state phase space is 
already occupied. The Pauli principle of the final state reduces the free cross-section to 
effective levels. Here we take different forms of NN cross-sections [cr] to understand the 
influence on fragment production. We shall here use a energy dependent cross-section due 
to Cugnon and also different constant cross-sections ranging between 20 and 55 mb [5,7]. 
The nature and strength of cross-section a is depicted as superscript to a. 



3. Results and discussion 



We here simulate the symmetric reactions of -;q,Au^' + 79 Au^' at incident energies of 
100, 400, 600 and 1000 MeV/nucleon and over complete range of the impact parameter. 
The system size and asymmetry effects will be analyzed by further studying the reactions 

ofioTVe^O + 13^^27, ,3^1^40 + 2l5c45, 36Xr84 ^ ^.TVfeSS and 54Xe"l ^ ^^^^139 

at energies between 20 and 150 MeV/nucleon. The clusterisation at first instance is made 
using the minimum spanning tree (MST) method, in which, nucleons are bound if Rcius ~ 
\n -rll <4fm. [1,5]. 
In fig.l, we display the multiplicity of intermediate mass fragments (IMF's) as a function 
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Figure 1. Comparison of average multiplicity of intermediate mass fragments (IMF's) 
with ALADIN data at incident energies of 100, 400 MeV/nucleon (top panel) and 600, 
1000 MeV/nucleon (bottom panel) as a function of impact parameter. The results are 
displayed using soft momentum dependent (SMD) interactions. 
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of the impact parameter for the reaction of igAu^^"^ + ygAw^^^ at incident energies 100, 
400, 600 and 1000 MeV/nucleon. Here soft momentum dependent interaction (SMD) is 
used with different cross-sections. From the figure, the multiplicity of IMF's is maximal 
at 100 MeV/nucleon for smaller impact parameters, which decreases with the increase in 
the impact parameter. On the other hand, one sees a rise and fall in the multiplicity of the 
IMF's at higher incident energies. The dynamics at 100 MeV/nucleon, is mainly governed 
by the mean field or by the density of the reaction as compared to other higher incident 
beam energies under consideration (e.g. 400, 600 and 1000 MeV/nucleon). The incident 
energy of 30 MeV/nucleon is the lowest limit for any semi-classical model, where the 
effect of Pauli-blocking is w 90%. Below this incident energy, quantum effects as well as 
pauli- blocking need to be redefined. There are no visible effects of different cross-sections 
at 100 MeV/nucleon. Due to the low excitation energy, central collisions generate better 
repulsion and break the colliding nuclei into IMF's, whereas for the peripheral collisions, 
the size of the fragment is close to the size of the reacting nuclei, therefore, one sees a very 
few IMF's. In contrary, a rise and fall can be seen at other higher incident energies. For 
the central collisions, the frequent NN collisions occurring at these energies do not allow 
any IMF's production, whereas, at peripheral collisions the energy transfer is from the 
participating matter to spectator matter is minimum, therefore, very few IMF's are seen. 
In all the cases, some effects of different NN cross-sections are visible at higher incident 
energies. The use of the momentum dependent interaction yields better comparison with 
ALADIN setup [18] for cr = 55 mb. This finding is in agreement with the results reported in 
ref. [1], where it was found that NN cross-section has sizable effect on reaction dynamics. 
Note that in these studies, static equation of state was used. 

Let us now analyze the above effects in asymmetric reactions . In fig. 2, asymmetric re- 
actions of loiVe^o + 13^^27^ ^g^^40 ^ 2^5-245^ aeKr^'^ + uNb^^ and54Xei3i + .^La^s 
are displayed as a function of beam energy at scaled impact parameter 6 = 0.3 (semi central 
collisions) using different cross-sections. Final results are also compared with the NSCL 
experimental data [19]. Due to no access to filters, no direct comparison with data could 
be made. These comparisons are just indicating the trend within theoretical framework. 
From the figure, it is clear that the trends of our calculations with soft momentum depen- 
dent (SMD) interactions are in good agreement with experimental data. Again different 
NN cross-sections fail to make any significant impact on the cluster dynamics. 
One further see that the multiplicity of IMFs in loiVe^" +13 AP''^ decreases with increase 
in the beam energies. For the reaction of ig^r^^ +21 Sc^^, similar trends emerge above 
55 MeV/nucleon, as in case of loA^e^^ +13 AP'^ . For energies below 55 MeV/nucleon, we 
see dominated role of mean field and hence increase in the intermediate mass fragments. 
For the rest of the reactions, namely, sgXr*^ +41 Nb^^ and 54X3^^^ +57 La^^^, the mul- 
tiplicity of IMF's increases with the increase in the beam energy and cross-section. One 
should note that in the first two reactions, incident energy is much higher compared to the 
last two reactions. 

Let us now understand how other technical parameters affect the fragmentation. In fig. 
3, we display in the upper part, the effect of width of Gaussian wave packet. We display 
the results with narrow width (L = 4.33 frn^) and broad one (L = 8.66 fm?)- We see that 
the variation of the width has sizable effect on clusterisation. A broader Gaussian results 
in extended interaction radius, therefore, binding more nucleons into a fragment. As a re- 
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Figure 2. The average IMF's multiplicity versus beam energy in the reactions of 

loiVe^" + x^Al^\ i8^r« + 2iSc«, 36^r«^ + 4iiV&^' and 54X6^2^ ^ ^^^^139^ 
The symbols represent the NSCL experimental results, while lines are representing the 
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Figure 3. A Comparison of multiplicity of IMF's by employing different model ingre- 
dients like Gaussian width (upper panel), clusterisation cut off distance Rcius (middle) 
and equation of state (bottom). These results are displayed using soft (upper, middle, 
bottom) and hard (bottom) equation of state with a'''' . 
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suit, the fragment turns much heavier than the upper limit of IMF's (i.e 65 for (Au+Au) 
reaction) and hence there is net reduction in the production of IMF's. This result is in 
agreement with the findings in ref. [1]. One should also note that larger Gaussian width 
will also result in more attractive nuclear flow and hence will push the energy of vanishing 
flow towards higher incident energies. In the middle panel, we display the effect of cut 
off distance Rcius on clusterisation. We display the results with Rdus = 4 fm and a 
narrower one 2.5 fm. We see that when we choose a narrow cut off distance, lesser number 
of fragments are formed. In the bottom panel, we display the role of different equations of 
state by simulating the reactions with hard and soft equations of state. We see that different 
equations of state do not alter the results. 

To further analyze the results, we display in fig. 4, the maximal number of intermediate 
mass fragments (Nimf) as well as corresponding impact parameter as a function of the 
incident energy for the reaction of rgAu^'^'^ + jgAu^^"^ . Here apart from the different 
Gaussian widths and clusterisation distance, we also display the results with different clus- 
terisation algorithms. Here we also use Simulated Annealing Clusterisation Algorithm 
(SACA) [l]to clusterise the phase space. In SACA approach [1], a group of nucleons 
can form a fragment if the total fragment energy/nucleon ( is below a minimum binding 
energy: 



N' 



c-E 



=1 






i¥'i 



< Lbe X Nf, (13) 

with Lbe = -4.0 MeV if N^ > 3 and Lbe = otherwise. In this equation, N^ is the number 
of nucleons in a fragment, P^ is the center-of-mass momentum of the fragment. The 
requirement of a minimum binding energy excludes the loosly bound fragments which 
will decay after a while. To find the most bound configuration among the huge number 
of possible fragmentation patterns, we proceed as follow: we start from a random config- 
uration which is chosen by dividing the whole system into few fragments. The energy of 
the individual clusters is calculated by summing over all nucleons present in that cluster 
using eqn.l3. Note that as we neglect the interaction between fragments, the total energy 
calculated in this way differs from the total energy of the system. If the difference between 
the old and new energy is negative, the new configuration is accepted. If not, the new 
configuration may nevertheless be accepted with a probability of exp(-AE/c), where c is 
the control parameter. This procedure is known as the Metropolis algorithm. The control 
parameter is decreased in small steps. This algorithm will yield eventually the most bound 
configuration. The identification of fragments using SACA is very time consuming. The 
algorithm searches the most bound configuration out of several million of different pos- 
sibilities. The major outcome of this method is that we can recognize the fragments at 
relative high density. 

From the figure, we see that different cross-sections yield quite similar trend. The ab- 
solute value varies with NN cross-sections as well as with other model ingredients. For 
example, we see that at 600 MeV/nucleon, the SACA method with 20 mb yields results 
close to the MST method with 55 mb. A broader Gaussian scales down the number and 
so is the case with reduced cut off distance. From this discussion, it is clear that the effect 
of larger cross-section (i.e. between 20 and 55 mb) is as large as effect of other model 
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ingredients. Even a change of the clusterisation algorithm affect the outcome sizeably. 
From this analysis, one conclude that it may not be possible to pin down the magnitude of 
cross-section from multifragmentation since other technical parameters such as width of 
the Gaussian, clusterisation range or even the change of the clusterisation algorithm alters 
the results in similar fashion when semi-classical model such as quantum molecular dy- 
namics is used. 



4. Summary 

In the present study, we focused on the comparative study of different model ingredients 
with semi-classical model, namely, quantum molecular dynamics (QMD) model. For this 
study, a hunt was made to compare our theoretical results of different cross-sections with 
experimental data. We also analyzed the results with model ingredients such as width of 
the Gaussian, clusterisation range and different clusterisation algorithms. We found that 
the effect of different cross-sections is of the order of the one obtained from the model 
ingredients. All model ingredients affect the fragmentation pattern in a similar fashion. 
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